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Abstract 

The isoscalar giant monopole resonance (ISGMR) and isovector giant dipole resonance 
(IVGDR) in finite nuclei are studied in the framework of a relativistic transport approach. The 
kinetic equations are derived within an effective nucleon-meson field theory in the Relativistic Mean 
Field (RMF) scheme, even extended to density dependent vertices. Small amplitude oscillations 
are analysed using the Relativistic Vlasov (RV) approach, i.e. neglecting nucleon collision terms. 
The time evolution of the isoscalar monopole moment and isovector dipole moment and the cor- 
responding Fourier power spectra are discussed. The excitation energies of ISGMR and IVGDR 
are obtained for spherical nuclei with various sets of Lagrangian parameters. 

In the case of 208 Pb we study in detail the dependence of the monopole response on the effective 
mass and symmetry energy at saturation given by the used covariant effective interaction. We 
show that a reduced m* and a larger 04 can compensate the effect on the ISGMR energy centroid 
of a much larger compressibility modulus K nm . This result is important in order to overcome the 
conflicting determination of the nuclear compressibility between non-relativistic and relativistic 
effective interactions. 

For the symmetry energy dynamical effects, we carefully analyze the influence of the inclusion of 
an effective isovector scalar channel, (5-meson field, with constant and density dependent couplings. 
We show that the <5-meson contribution, keeping fixed the equilibrium 04 value, leads to a change 
in the isoscalar and isovector response of finite nuclei which results in a reduction in the centroid 
energy of ISGMR and IVGDR for 208 Pb. All that in fact reveals the relevance of the slope (or 
pressure) of the symmetry energy at saturation on the ISGMR and IVGDR modes for neutron- 
rich systems. 

Density dependent vertices are not much affecting our conclusions. Following as a guidance some 
extended dispersion relations in nuclear matter, we see two main reasons for that, the smoothness 
of the density dependences around saturation and the presence of compensation effects coming 
from rearrangement terms. 

PACS numbers: 21.30.Fe, 21.65. +f, 24.30. Cz, 21.10.Rc 

Keywords: Relativistic Transport, Collective Response, Nuclear Compressibility, Symmetry Energy, Effec- 
tive Hadron Lagrangians. 
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I. INTRODUCTION 



In this work we study the isoscalar giant monopole and isovector giant dipole oscillations 
in 208 Pb using the Relativistic Vlasov method (RV) in a Relativistic Mean Field (RMF) 
theory with constant and density dependent meson-nucleon vertices. We analyze various 
RMF models with different coupling parameters which lead to a different nucleon effec- 
tive mass (m*), compressibility modulus (K nm ), symmetry energy (04), keeping fixed other 
nuclear matter saturation properties. 

Among those models we focus on the ones where the scalar isovector channel has been 
introduced by the coupling to an effective 5-meson field jlj. It has been shown that such 
inclusion has important effects on the equation of state (EOS) and phase diagram of asym- 
metric nuclear matter (ANM) jji . as well as on the reaction dynamics with exotic nuclei, 
see the review [2jj]. In fact the 5- meson brings contributions to the slope and curvature of 
the symmetry energy and to the neutron-proton effective mass splitting. In particular the 
influence of this coupling on the collective response of (ANM) appears important, as shown 
in a linear response approach in ref. Q]. For the comparison of the results we have also con- 
sidered the well known NL3 0, Q| parametrization, very successful for finite nuclei structure 
calculations, and some density dependent vertex models (f| Q| beyond RMF. In this way 
we aim to pin down the sensitivity of isoscalar and isovector collective oscillations on m*, 
K nm and a 4 , focussing on the interplay between isoscalar-isovector couplings. In particular 
we will see that the covariant structure of the isovector part of the effective interaction will 
clearly show up even in the isoscalar monopole mode in 208 Pb. 

Relativistic mean field models, with constant jg, [| and density dependent [n]] meson- 
nucleon couplings, have been applied to the description of collective excitations of atomic 
nuclei in the framework of Time-Dependent Relativistic Mean Field (TDRMF) and in the 
self-consistent relativistic random phase approximation (RRPA). The monopole predictions 
of nonrelativistic Hartree-Fock plus random phase approximation (RPA) calculations, using 
both Skyrme and Gogny effective interactions, seem to indicate that the value of K nm should 
be in the range 210-220 MeV [ill . 12 1. In relativistic mean-field models on the other hand, 



results of both time- dependent and RPA calculations suggest that empirical GMR energies 
are best reproduced by an effective force with K„. m 250-270 MeV This difference has 

been pointed out as a serious open problem [lOj. A possible solution to this ambiguity was 
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suggested in Ref. [3] as an interplay between the effective mass and compressibility within 
a discussion of isoscalar response of nuclear matter. Here we show the validity of this 
interpretation, joint to a density dependence of the symmetry energy , also for the spherical 
finite n-rich nuclei, like 208 Pb. Recent studies are in fact pointing to the same isospin effect, 



see refs 



. |l5 
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17| . In this respect we have even analyzed the influence of the inclusion of 
the 5-meson on both isoscalar and isovector response in spherical n-rich nuclei. 

The introduction of density dependent vertices, following the Dirac — Brueckner — 
Hatree — Fock estimations, is not showing relevant new effects. A nice analitical inter- 
pretation of this result, based on a linear response theory for nuclear matter, is finally 
presented. 

II. OUTLINE OF THE MODEL 

The dynamics of collective vibrations in spherical nuclei is studied in the framework of 
the relativistic Vlasov (-RV") transport equation, which describes the dynamical evolution of 
a semi-classical phase space distribution function f(x,p) under the influence of the nuclear 
mean field. Binary collisions between nucleons will be not considered here. Thorough 
derivations of the RV transpor t eq uations from an effective hadron-meson field theory, 



can be found elsewhere 
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2l\ . The RV equation reads (i = p,n): 
\pT% + + m^m^V^p*) = (1) 

whith the field tensor 

FT = - <9"Ef (2) 

and effective masses and kinetic momenta m* and p* 11 , respectively, specified below. The 
particles obey the mass-shell condition 

P*% - mf = 0. (3) 

Thus, from the temporal knowledge of the phase space distribution function one can 
calculate the time evolution of physical quantities such as densities and fields. We remind 
the meaning of the the Wigner matrix (in the Spinor space) 
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F a p(x,p) = <F a p(x,p)> 

= j^r I d ' R e ~ iP ^ < + l R)Mx ~\ R)> - (4) 

^From the above definition it follows that single-particle operators can be expressed as 
(spinor indices will be omitted for simplicity) 

<0>= J d A x J d A p tr (dF&pfj (5) 

where the trace runs over spin and isospin indices. The scalar density and the vector current, 
for example, assume the form 



p s {x) = < ipip >= J d 4 p tr(F(x,p)) (6) 
j^x) = < ipj^ip >= / d 4 p tr(7 M F(x,p)), (7) 



and will be used to calculate the different Lorentz components of the mean field potential. 

The nuclear mean field U is characterized in a RMF theory by means of self energies in 
the form U = T, s — 7 M S M • • • (higher contributions are usualy neglected due to symmetry 
properties of nuclear matter). In a Non-Linear (QHD — NL) model with isoscalar scalar 
and vector meson fields a and u and with the inclusion of the isovector channel through the 
exchange of the virtual charged ^(scalar) and p(vector) mesons, the mean field approximation 
leads to self energies which are related to the expectation values of the combination of 
isoscalar and isovector fields with coupling constants g a , g w , g p and gg. The scalar and 
vector components of the self energies are generally given by 

. + proton (i = p) 
g^(x)±g p b»(x) { P V FJ (8) 

— neutron {i = n) 

f proton (i = p) 
g a a(x) ± g s 5(x) { (9) 

neutron {i = n) 

with the expectation values of the fields self-consistently calculated, see later Eqs. (fT??|) . 

The self energies characterize the in-medium properties of the nucleons inside the hadronic 
environment in terms of kinetic momenta and effective masses 

vT = v i t-n (io) 
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m* = m - S si (11) 

The density dependence of the mean field, i.e. the density behavior of the self energies, 
depends on the coupling choices of the RMF model. Here we consider different parametriza- 
tions within the Non-Linear (QHD — NL) effective field approach and even extended to the 
Density Dependent Hadronic (DDH) mean field theory jy, Q, 22]. In table H] the values 
for the different coupling constants and the non-linear parameters for different sets of Non- 
Linear Walecka (QHD — NL) models are presented, for details see refs. 
Their corresponding nuclear matter saturation properties are given in table |H] In the DDH 
models the density dependence of the coupling constants is chosen in order to reproduce 
microscopic Dirac — Brueckner — Hartree — Fock results beyond the RMF picture, see 



the discussion in ref. 



221. 



The choice of models with rather different nuclear matter properties has been done on 
purpose, in order to investigate the role of the effective masses, compressibility and symmetry 
energy on isoscalar (monopole) and isovector (dipole) oscillations in neutron-rich nuclei. We 
have compared various Non-Linear RMF parametrizations, in particular the Giessen sets 
NL1 — G, NL2 — G, [^J, extended also to simulations of relativistic heavy ion collisions, and 
the widely used NL3 set , successfully applied in finite nuclei studies. The different 
treatment of the iso-vector part of the mean field (competition effects of the repulsive p 
field and the attractive 5 meson) is analysed in detail using the (NLp, NLp5) parameter 
sets pj. The same analysis is performed within density dependent coupling models, the 
parametrizations DDH3p and DDH3pS of ref. (3] have been used. 

The dynamics of collective vibrations in spherical nuclei is simulated in the framework 
of the relativistic Vlasov equation (JIJ. Its numerical realization is based on the standard 
test particle method, where the phase space distribution function f(x,p) is represented by 
a finite number of test particles of a covariant Gaussian form (Relativistic Landau- Vlasov 
(RLV) method). A detailed description of the RLV method is given in Ref. j^jj], here this 
model will be briefly discussed. The use of a Gaussian shape for the test particles has the 
advantage of smooth distribution functions, but maintaining an accurate determination of 
local quantities, particularly important near the nuclear surface. 

The covariant Gaussians in the four- dimensional Minkowski space are defined as 
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TABLE I: Coupling parameters in terms of f\ = {gi/rrii) 2 for i = a, to, /, = (gi/2nii) 2 for p,5, 
A = ajg^ and B = b/g^ for various Non-Linear RMF models using the p and both p and 5 mesons 
for the characterization of the isovector mean field. In the DDH models the coupling functions 
are explicitly density dependent 
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Parameter 


NL1-G 


NL2-G 


NL3 


NLp 


NLp5 


Ufm 2 ) 


6.146 


9.300 


15.739 


10.330 


10.330 


fUfm 2 ) 


3.611 


3.611 


10.530 


5.420 


5 420 


fp(f™ 2 ) 


1.200 


1.220 


1.339 


0.950 


3.150 


fsifm 2 ) 


0.0 


0.0 


0.0 


0.0 


2.500 


A{fm- 1 ) 


-0.123 


0.0824 


-0.01 


0.033 


0.033 


B 


0.169 


-0.0012 


-0.003 


-0.0048 


-0.0048 


TABLE II: Nuclear matter saturation properties in the different RMF 


models. 


Property 


NL1-G 


NL2-G 


NL3 


NLp, p5 


DDHp, p5 


E/A(MeV) 


-16.0 


-16.0 


-16.3 


-16.0 


-16.0 


Po(fm- 3 ) 


0.145 


0.145 


0.148 


0.160 


0.153 


m* jm 


0.83 


0.83 


0.60 


0.75 


0.55 


K nm (MeV) 


380 


210 


272 


240 


240 


a 4 {MeV) 


30.62 


30.62 


37.40 


30.50 


33.40 



G(x;£i) = / drg(x-Xi(r)) 

J — oo 

/ — CO 
dr exp((a; - x^t)) 2 /w 2 )5[(x^ - x^(t))u^{t)} (12) 
-oo 

where £j denotes the world line of the particle j as a whole, r refers to the eigentime of the 
test particle and a is the normalization constant. In the four- dimensional momentum space 
a gaussian weight of a test particle is defined by 

g(p* - p*(r)) = a p exp((p* - p*(r)) 2 /w 2 p )6[ P ; P r(r) - mf] (13) 

where the center of the gaussian is assumed to be on-shell i.e. p* ifjL = m*Ui^, uf = 1, whereas 
the free momentum is generally off-shell. The effective mass of the particle is taken as 
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m* = m*(xi(r)). The norm of the gaussian is calculated in the rest frame of the particle to 
be m*' 1 with a p = (y/7rw p )~ 3 , where w and w p are the test particle widths in position and 
momentum space. 

With the Gaussians of (jl2piHj) the phase space distribution function f(x,p*) is expressed 

as 

■y A-N f+oo 

ffap*) = / drg(x - Xi(r))g(p* -p*(r)) 

„•_ 1 J— oo 



i=l 
1 



A-N 



N(ttww p ) 3 



1=1 



S y^ j j dr exp((x — Xi(r)) 2 /w 2 ^ 



x exp((p*-p*(r))7^) 



(14) 



where N is the number of test particles per nucleon. Scalar densities p s and baryon currents 
jn follow from the phase space distribution 

A-N 

N 

A-N 



^■ r) = / ~~ dT J^} X ,LM x-Xi(T)) (15) 

i=l * / -° 



1 E 



i=l 



m*(xj(r))' 
dr^x - Xi(r))M iM (r). 



The equations of motion for the test particle trajectories are given by 



(16) 



d_ 
dr~ 
d_ 

dr 



u?(t) 



ut(r) 



A-N „ o 

i 

m*(xi) ^— ' w 2 Vm 
da(x 



I \ 5 v 2 [ 2 U iv(R-'j 

\Xi) ^— ' ttr Lm^ J 



dp s 



exp(R 2 (xi)/w 2 ) 



± 



2 ^ r ^ 



Am 2 



-j exp(.R|(x i )/w 2 ) 
N(y/irw) 3 



A-iV 



m*(xi) w 2 j= f^ +l Umj 



] exp( J R|(x i )/w 2 ) 



(17) 



(18) 
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with 



R?(x) ee - <(r)) - (x„ - MrJKWM, 



projection of the vector (x — Xi(r)) on the hyperplane perpendicular to Ui(r), 24]. 

Here the equations for velocities, rather than for momenta, are given, within the assump- 
tion that the particle accelerations are small. In this equation particles are propagated by 
their respective eigentimes and so corresponding time coordinates x^t) can differ. In order 
to solve the pro blem of different time coordinates a system time has been adopted for the 



propagation 



241- 



The scalar and vector meson fields, determined by the scalar density p s and the baryonic 
current j^, respectively, result from the solution of the corresponding equations in the Local 
Density Approximation (LDA) 



mla(x) + Ba 2 (x)+Ca%x) = g aPs (x) = g a J d^*'-^f(x t p*) 

= ^j„(x) = / d 4 p*p* f{x,p*) 
ml m^ J 

W = oi^) = A / d*p%f 3 (x,p m ) 
p p J 

with f 3 = f p - f n . 

An important issue for the description of low energy excitations within the RV equation 
is an appropriate initialization of finite nuclei, before starting the phase space evolution of 
the distribution function with Eq. (|18|). The ground state of a spherical nucleus is randomly 
initialized by means of the test particles with the covariant Gaussian shape in position and 
momentum space. After the first step of a randomly obtained initial distribution, a fit 
procedure in coordinate space is performed according to realistic Thomas-Fermi reference 
density distributions. During this procedure the proton and neutron distributions are fitted 
separately by satisfying the empirical values for the asymmetry and the surface thickness 
parameters, respectively. In the Thomas-Fermi calculation the scalar densities are deter- 
mined by solving self-consistently the equations for the effective masses m* n . With this 
initialization of a given spherical nucleus the temporal evolution is described by the RLV 
equation (fTH|) . We have used 100 test particles per nucleon for the transport descriptions, 
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which yields a smooth distribution function with a very good energy conservation. 

For giant resonances the time dependence of collective dynamical quantities is not peri- 
odic, since giant resonances are generally not stationary states of the mean-field Hamiltonian 
0| . For non-relativistic models it has been proven that the results from the Vlasov equation 
are identical to results from full quantum mechanical calculations (TDHF) 25[. Since in 
the non-relativistic frame RPA is the small amplitude limit of a TDHF calculation, one can 
expect that the results from the relativistic Vlasov equation are comparable to those from 
relativistic RPA calculations of the type discussed in Ref. Q]. In the small amplitude limit 
the energy obtained from the frequency of the oscillation coincides with the experimental 
energy of the collective oscillation. 

The collective dynamical variables that characterize nuclear vibrations are defined as 
expectation values of single particle operators in the phase space representation. For the 
isoscalar monopole vibrations, the time dependent monopole moment is defined as 

<r 2 (t)> = j J d 3 xr 2 Ux). (20) 

In this work we have applied the RLV method to isoscalar monopole and isovector dipole 
oscillations in 208 Pb. The excitation of an iso-scalar monopole oscillation of the initialized 
nucleus in its ground state is modeled by a radial expansion. This is done by introducing a 
new coordinate 

r m on = (1 + a)r (21) 

for each test particle. Here a is a scaling parameter, a = 0.1 fm has been used, and the 
deformation of protons and neutrons is in phase. 

For the isovector dipole oscillation the following operator has been applied [2"i 

N Z N 

«T = nTz £ r ^ - nTz £ (22) 

p=l n=l 

This means an out of phase shift along the z-axis between protons and neutrons. We have 
used a scaling parameter of 1 fm according to Eq. ()22|) . which causes a center of mass 
separation between protons and neutrons, while keeping unchanged the center of mass of 
the whole system. 
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FIG. 1: (Left) Time-dependent isoscalar (r 2 ^ monopole moment and (right) the corresponding 
Fourier power spectra for 208 Pb, as obtained from transport calculations with different choices of 



the nuclear mean field: (NL1 — G, NL2 — G, NL3) non-linear parametrizations of Refs 
respectively. The corresponding excitation energies of the ISGMR are indicated in the right 
pannels. 



III. ISOSCALAR MONOPOLE OSCILLATIONS 



The study of Isoscalar Giant Monopole Resonances, ISGMR, in nuclei provides an impor- 
tant source of information on the nuclear matter compressibility. The complete experimental 
data set on isoscalar monopole resonances has been analyzed by Shlomo and Youngblood 
j^J. In fact, within a semi-empirical macroscopic approach of the systematics of GMR, it 
turns out that starting from finite nuclei data the nuclear matter compression modulus K nm 
can be only fixed within the rather wide range 200 — 350 MeV. A systematic theoretical 
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work, based on non-relativistic RPA methods with realistic effective interactions |l2j ]. was 
supporting the lower region of the above range, around 210 — 220 MeV. At variance rela- 
tivistic approaches, within the same small amplitude limits, were pointing to some definitely 



larger values, 250 — 270 MeV lOj. We address the problem studying in a relativistic frame 
the ISGMR in 208 Pb, that shows a well established Giant Monopole Resonance (GMR) 
at 13.7 ± 0.3 MeV. Moreover 208 Pb is interesting since it represents a neutron-rich nucleus, 
where we can even reveal isovector effects. 

As discussed in the previous section, the iso-scalar monopole oscillations are analyzed 
within the relativistic Vlasov transport equation. The nuclear mean field is evaluated in the 
Non-Linear versions of the Walecka model and in the DDH approach, see Tables and [nj 
We will discuss the dependence of the GMR on the compression modulus, e.g. by comparing 
the Giessen parametrizations (NL1-G and NL2-G) with a similar value for m* but different 
K nm , and on the effective mass (for fixed K nm ), e.g. by comparing the parametrizations of 
Liu et al. with the DDH model. The latter will allow also some comments on the effect 
of a density dependence of the couplings around the saturation point. These analyses are 
comparable with similar studies in the framework of a time dependent relativistic mean field 
model, see Refs. [J HQ, and within non-relativistic RPA calculations of the nuclear 
matter compressibility using Gogny effective interactions, see Ref. |12|. 

The influence of a different treatment of the iso-vector part of the mean field will be 
particularly discussed either paying attention to the differences of the various interactions, 
see Table ITU or more specifically focussing on the effects of the iso-vector, scalar 5 field, i.e. 
by comparing the results of parametrizations, like NLp, NLpS and DDHp, DDHpS, with 
and without the introduction of the 6 meson. We note again that the inclusion of the S meson 
in the interaction leads to a stiffer symmetry energy (around saturation) and to a splitting 
between the (Dirac) effective masses of neutrons and protons, see a review in ref. \2]. Also 
the vector component of the iso-vector self energy is modified due to the enhancement of 
the vector, iso-vector p meson field roughly by a factor of 3. 

Fig. 1 shows the dynamical evolution of the monopole moment of the excited 208 Pb nu- 
cleus for those models with fixed effective mass (NL1—G, NL2—G) and different compression 
modulus. The ISGMR oscillation is modeled by Eq. ()20|). The numerical simulation of 
the RV equation is performed with the test particle method of the previous section. The 
limited number of test particles per nucleon (100 in our case) leads to numerical fluctuations, 
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which can be seen in Fig. ^ in terms of spurious oscillations with very low frequencies (first 
peaks in the Fourier energy spectra) and in terms of a partially non-periodic evolution of the 
monopole moment. This is a general non-trivial feature of transport de scrip tions with test 
particle methods which leads to an artificial damping of the excitation 30]. However, due 
to a rather good stability of the nucleus (until several 100 fm/c) and a good separation of 
the ambigous excitation energy of the numerical oscillations, we conclude on the reliability 
of the results presented in the following. 

The excitation energy of the ISGMR oscillation is very sensitive to the nuclear matter 
compressibility, which is a well known fact. In particular, the excitation energy increases 
with rising compression modulus (from NL2 — G to NL1 — G)), as expected. 

The experimental value of 13.7 ± 0.3MeV can be reproduced satisfactory by the NL3 
model (13.39 (hu)). This result is consistent with the analysis of Ref. 4J], again as an 
important check of the transport calculations. 

We notice however that the NL2 — G model, with K nm = 210 MeV, leads to a centroid 
excitation energy (11.78 (hoo)) not too far from that of the NL3 model, which has a much 
larger compressibility, K nm = 272 MeV. Considering the other saturation properties of 
the NL3 force, see Table |HJ this fact is pointing to two interesting implications: i) The 
dependence of the monopole frequency on the effective nucleon mass; ii) The dependence on 
the symmetry energy, the 04 parameter, that one could expect since the 208 Pb is a neutron- 
rich system. We will discuss in detail the two effects in the following. We note that, as a 
straightforward consequence, the claimed discrepancy on the nuclear matter compressibility 
modulus K nm between relativistic and non-relativistic models, see [lflj . can be eventually 
overcome. 

1. Monopole frequency and effective mass 

A linear response theory for nuclear matter within the RMF frame can give some inter- 
esting hints. We can derive a dispersion relation, see refs. 



E* F 
1+ ' 



3 k'jp 



-9f»^(l-U^ p s ) Pb] <p(s) = . (23) 
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where s is the dispersion parameter s = ^ = and </?(s) is the usual Lindhard function 
of the Landau Fermi Liquid theory: 



s + 1 



+ ^7TS0(1 - S) 



8-1 

Here the is the potential part of the nuclear matter compressibility. 

From Eg. (|23jl we see that the "restoring" force for monopole oscillations is given by an 
effective compressibility which is reduced for larger values of the u meson coupling constant. 
However / w can assume very different values depending on the chosen value for effective 
masses m*. This is easy to understand since in the RMF limit the saturation binding 
energy has the simple form 

E/A(p ) = Ep + f^po - m N 

where is the bare nucleon mass. Thus in order to have the same saturation values of po, 
E/A(po), when we decrease m* we have to increase f w . We then come to the natural conclu- 
sion that two EOS with different effective masses at saturation, even if the compressibilities 
are the same, are expected to have different dynamical monopole response. In the NLG — 2 
vs. NL3 comparison discussed here we clearly see the interplay between compressibility and 
effective mass: we can get similar monopole energies increasing the compressibility while 
decreasing the effective mass, see Table HU 

This appears to be a quite general feature, present also in non-relativistic approaches, 



12j where the RPA systematics of the Gogny forces is shown: the 



see the Fig. 7 of the ref. 
208 Pb breathing mode energy is not much changing if at the same time we increase the NM 
compressibility and decrease the effective mass. A similar trend has been suggested in a 
recent work on the nuclear compressibility within the nonrelativistic frame [itI ]. 



2. Monopole frequency and symmetry energy 



It is well known that the equilibrium properties of nuclear matter are changing with 
isospin asymmetry, in particular the saturation density and the corresponding EoS curvature 
Sect. 2 and refs. therein. For the compressibility shift we have, after some algebra: 



, see 



AK nm (I) = 9p c 



d 2 2 d 
^° dp 2 dp. 



-syvn 



(P) 



P=PO 



[K sym -6L]I 2 < 0, 



(24) 
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FIG. 2: (Left and right pannels) Time-dependent isoscalar (r 2 ^ monopole moment and (middle 
pannels) the corresponding Fourier power spectra for 208 Pb. Transport calculations with the pa- 
rameter sets are shown: (panels on the top) Liu — RMF with (p) dashed line (right) and with 
(p + 8) solid line (left), (panels on bottom) Density Dependent model (DDH) (lower panel) with 
(p) dashed line (right) and with {p + 5) solid line (left). The excitation energy of ISGMR for each 
parameter set is explicitly indicated in the Fourier power spectra and for comparison the results 
with (p + S) and (p) for both parameter set are shown together. 

with / = (p n — p p )/p, asymmetry parameter. We note the interplay between slope, L = 
3pot' S ym(p), and curvature, K sym = 9poe" ym (p) of the symmetry energy at saturation. The 
asymmetric matter becomes softer since the shift is in general ne gati ve, due to the dominance 
of the slope term L 0, see Q and the recent discussion in [22|. Thus, in 208 Pb, n-rich 
system, / = 0.21, asymmetry can affect the the isoscalar monopole oscillations, as also noted 
in refs. ^3, 16, 3| within relativistic and non-relativistic frames. 
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In our NL3 vs. NL2 — G comparison the difference on the values, see Table |Hj auto- 
matically implies a difference on the slope parameter L, since in both models the potential 
symmetry energy is coming from the effective p-meson coupling which leads to a linear p#- 
dependence. So the larger 04 of NL3 means a larger slope at po an d so a larger reduction 
of the compressibility in 208 Pb, see Eq. (J24|) . Our conclusion is that in NL3 both effects, 
smaller effective mass and larger 04, are almost compensating the much larger K nm value, 
finally leading to a monopole frequency in 208 Pb not much different from the NL2 — G one. 

We have continued the analysis of the symmetry contribution studying the effect of 
a scalar isovector channel with the inclusion of 5 meson both in nonlinear and density 
dependent models. The results are shown in FigJ21 The time history of isoscalar monopole 
moment and its Fourier spectrum are shown for the Liu sets in the upper panel with isovector 
p + 5 interaction (NLpS, solid line) and with p interaction (NLp, dashed line) and similarly 
with density dependent models (lower panel). For both interactions although the time 
evolution shows a very similar behaviour the power spectrum presents a net reduction in 
the energy centroid when one includes the 5 meson. This is a very nice indication of the 
dominance of the "Slope" of the symmetry energy on the compressibility shift. In fact when 
the 5-channel is included we have a clear increase of the symmetry energy slope L around 
saturation, see the discussion after Eq.(6-19) of ref. e.g. it results about 20% larger in 
the Liu parametrizations. 

Our discussion suggests that although the inclusion of S meson does not produce impor- 
tant effects on the ground state of finite nuclei 3], it has an interesting influence on the 
collective excitations of charge asymmetric finite nuclei, in particular for 208 Pb. This ap- 
pears a good suggestion for new experiments aiming to a better determination of the poorly 
known slope (and curvature) of the symmetry term around saturation. 

In general the DDH monopole frequencies are systematically below the corresponding Liu 
ones. However this appears mostly a joint effect of different effective masses and symmetry 
energies, as already discussed for the NL3 vs. NL2 — G comparison. Therefore the density 
dependence of the coupling constants seems not largely affecting the monopole response. In 
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fact this can be expected from the rather smooth behaviour around p , see Fig.l of ref, 
A more detailed study is presented in the Section V. 

Finally we like to note again that, although the Liu's NLp parametrization has a smaller 
K nm compared to NL3, it produces roughly the same monopole main frequency due to a 
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larger effective mass joint to a smaller 04. Of course we cannot compare the NLp8 results 
since the 5 channel is absent in NL3. 



IV. ISOVECTOR DIPOLE OSCILLATIONS 

In the literature only Relativistic RPA calculations have been performed so far to study 
the well known isovector giant dipole resonance It has been reported that 

for calculations of higher multipole modes, other than monopole, the response of a nucleus 
is difficult to evaluate in the time- dependent RMF [14|. The reason for this is that, since 
rotational invariance is broken and the differential equations have to be explicitly solved at 
each time-step on a two dimensional mesh in coordinate space, it becomes difficult to keep 
the solutions stable for the very long times that are necessary for a good accuracy. The 
problem is overcome in the present RBUU simulation, since it is possible to study the time 
history of IVGDR in a time dependent frame as already done in the ISGMR case. 

The experimental IVGDR energy in 208 Pb is well known as 13.5±0.2 MeV j^j]. Since in 
the isovector channel we mainly want to study the effect of the inclusion of isovector-scalar 
couplings, we only present results of the models that are parametrized with and without 5 
meson. In Fig. El we present the 208 Pb IVGDR oscillations and the corresponding Fourier 
transforms within models that include isovector-scalar channel namely Liu (upper panel) 
with p + 5, solid line(left) and p, dashed line (right) and similarly for Density Dependent 
parameter sets (lower panel). The Fourier transforms show a good single frequency dom- 
inance of the isovector mode. We observe that the DDH models, with 04 = 33.4 MeV, 
systematically give a larger resonance energy compare to Liu — RMF sets, with a\ = 30.5 
MeV. Moreover in both models a clear reduction of centroid energy is observed when the 
5-meson is included. 

It is a well known fact that the dynamics of IVGDR is rather sensitive to the symmetry 
energy of the corresponding model which is acting as a kind of restoring force parameter. The 
drawback in previous relativistic models is the one-to-one correspondence between 04 and 
IVGDR energy Q], and so it is difficult to discriminate among the different interactions. 
This is not the case in our analysis. The new important conclusion that can be drawn from 
our results is that the IVGMR dynamics is also sensitive to the more microscopic covariant 
structure of the symmetry term, i.e. to the interplay of various isovector channels. The 
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FIG. 3: Time-dependent isovector dipole moment (center of mass displacement of protons and 
neutrons (Rp C rn) — {R n cm)) and the corresponding Fourier power spectrum for 208 Pb. The param- 
eter set of the effective Lagrangian is "Liu" (upper panel) with (p) dashed line (right) and (p + 5) 
solid line (left), and Density Dependent model (DDH) (lower panel) with {p) dashed line (right) 
and (p + 5) solid line (left). The excitation energy of ISGMR for each parameter set is explicitly 
stated on the Fourier power spectrum and for comparison the results with (p + 5) and (p) for both 
parameter sets are shown together. 



physical interpretation of this result can be given in terms of the isovector response derived 
Ref.Q, where it is shown that the potential part of the symmetry energy explicitly appears 
in the dispersion relations with an correction term having a definite f p , fg structure: 
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(25) 



Note the similarity with the corresponding dispersion relation for isoscalar modes Eq.(|23J). 
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in particular the parallel role played by comressibility and symmetry energy in the two 
collective degrees of freedom. 

Now we can easily have interactions with the same value at normal density but with 
rather different isovector dipole response. In fact when including the 5 channel we have to 
increase the f p coupling in order to have the same symmetry energy at saturation (5, scalar 
field, is attractive in the isospin degree of freedom, see the discussion in Sect. 6 of ref.Q). 
Then the isovector dipole " effective restoring force" (coefficient of the Lindhard function in 
the Eq.(J2Hj) will be reduced. 

A similar effect has been pointed out in a detailed non-relativistic Skyrme — RPA study 
of the Giant Dipole Resonance in heavy nuclei ( 208 P6) using effective interactions with 
various isovector terms [3J|. A separate sensitivity of the average resonance frequencies on 
the symmetry energy 04 and on its slope has been found. In a covariant scheme we can see 
from Eq.(|25p that such behaviour can be achieved only by using two isovector fields, at the 
lowest order. This result shows more generally that a dynamical observable can be more 
sensitive to the microscopic structure of the isovector interaction than static properties. For 
instance in a careful study of the neutron distributions, [35], it is clearly shown that these 
"equilibrium" observables are almost equally correlated to value, slope and curvature of the 
symmetry term. 

Finally we see that even for the isovector dipole we cannot reveal specific contributions 
related to a density dependence of the coupling constants, in the sense that all the observed 
differences can be accounted for just in terms of symmetry energies and of isovector channels. 
We must say that we have not really performed an accurate analysis of this point, e.g. 
comparing with RMF models with exactly the same saturation properties. Enlightening 
analitical results can be derived within a nuclear matter linear response theory, as shown in 
the following Section. 

V. NUCLEAR MATTER RESPONSE WITH DENSITY DEPENDENT VER- 
TICES 

In previous Sections in order to explore the relevant physical quantities affecting the 
energy centroid of the collective modes, we have compared results from various QHD models 
and parametrizations. Furthermore to have a guidance on the observed effects in finite nuclei, 
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we took advantage of the plain dispersion relations for collective modes in nuclear matter 
studied in Ref. 3( in the contest of QHD — NL models. However in the DDH approach, 
where meson-nucleon couplings depend on the vector (baryon) density, the relation between 
the coupling functions and the compressibility can be expected to be modified together with 
the dispersion relations for isoscalar and isovector modes. 

By a comparison of results with the different parametrizations, we have argued that 
the density dependence of meson-nucleon couplings do not carry specific contributions at 
least to a large extent. In order to corroborate such a statement we briefly discuss the 
pertinent modifications of the thermodynamic quantities together with the linear response 
theory in nuclear matter, showing that in DDH models most of the effect is reabsorbed 
by the compressibility and large corrections may be expected only for very strong density 
dependence of the couplings. 

We remind that the DDH model has all the meson-nucleon couplings dependent on 
the isoscalar vector (baryon) density, defined as g = ^/J^J^, with j M defined as in Section 
ITT1 Such a density dependence leads to rearrangement contributions that affect the vector 
self-energy adding a term that at mean field level reads as: 



df u 2 d fv 2 , 9 f P ■ -v^dfs 2 \ 3ft 



^t'-W'^t^t&yz (26) 

The effect of rearrangement terms on the pressure can be find in refs. 0] • More relevant 
for our discussion is how they enter the compressibility: 

K = 3-£ + 9 [f uPB - U-et^Pb +9 -^p| - i^^-PsPB + 9 -^p B (27) 



Ep \ E* dps J \dps dps E* J dp B 

with 

d^o ( df u df a dp s \ 1 / d 2 f w 2 d 2 f a 2 \ 

dp B \dp B dp B dp B J 2 \dp B dp% J 

we can see that Eq. (j27|) reduces to the standard formula of QHD with costant couplings 
when the coupling are density indipendent, see Ref. Q|. 

The linear response in nuclear matter can be derived along the same lines of Ref. j^j , but 
one as to redefine the vector self-energy E M , defined in Eq.(JBJ), adding the rearrangement 
term and consequently the effective kinetic momenta, defined in Eq. dlUj) . Once also 
the equilibrium Wigner matrix is consistently redefined, on the same fashion of Ref. [3] a 
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dispersion relation for isoscalar modes in symmetric nuclear matter can be written as: 



3fc| [" nm " JW Ef y J " E$ ra ' E* F dp B 

dm 



(p{s) = Q. (29) 



if terms associated with J^r are discarded, or in other words, if the safe approximation 
jj^- — ^ is taken. We can see that of the derivative of the coupling function are reabsorbed 
in the K^. The difference between Eq.(J2HJ) and Eq. (j29j) is the last term. With some algebra, 
within the same approximation scheme, the contribution of density dependent terms can be 
reduced to a variation of the a-coupling, from f a to f a — P_b, in the correction term of 
Eq.(|2^j). Due to the other quenching factors and to the smoothness of the / ct (pb) function 
around po, we can expect an overall variation of a few percent in the dispersion relation. 

We note however that now the fuj{po) is not exactly the same of QHD — NL (with scalar 
non- linear terms), because the rearrangement terms affect also the relation for the binding 
energy, that in DDH models is given by: 

E/A(p ) = E* F + Upo + ^-m N (30) 

From Eq. (|2l)|) we see that is given by the difference in the density slope between the scalar 
and vector field. Therefore it may not be discarded in general, but is usually subdominant 
due to the similar density dependence of the scalar and vector coupling functions around 



saturation density, see Fig.l of ref. 22]. 



The expression for the symmetry energy is not modified in DDH models respect to 
QHD — NL [1], because the couplings depend only on the isoscalar vector density. Also the 
dispersion relation for the isovector mode in symmetric nuclear matter is not affected by the 
density dependence of the couplings and the relation written in Eq. (|23j) is still valid. However 
in asymmetric nuclear matter new terms appear, that are proportional to the asymmetry 
/ and to the derivative of the p— like and 5— like couplings. Therefore in DDH models for 
exotic nuclei there is the possibility to have a modified relation between symmetry energy 
and dipole excitation, however an exhaustive study of such effects goes beyond the scope of 
this paper. 
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VI. CONCLUSIONS 



Aim of this work has been to study effects on collective excitation properties of asymmet- 
ric spherical nuclei of the detailed structure of the nuclear EoS in the isoscalar and isovector 
sector. Particular attention has been put on the dynamical implications of the inclusion of 
a scalar isovector contribution due to an effective 5 meson coupling. 

We have described the dynamics of the isoscalar monopole and isovector dipole oscillations 
in a relativistic transport approach based on a nucleon-meson effective field interaction. 
We have evaluated the time-evolution of the oscillations obtaining then the corresponding 
excitation energies from the power spectra of the modes. The applications are to the n-rich 
nucleus 208 Pb, using a suitable choice of different effective Lagrangians. 

For the monopole mode we show an interesting m*, effective mass, dependence of the 
centroid energy. This effect, joint to a symmetry energy contribution, can account for the 
claimed ambiguity on the difference in compressibility moduli, that produce the experimental 
energy of ISGMR in 20S Pb , between relativistic and nonrelativistic microscopic models. 
From the influence of the 5 meson, just a systematic reduction of the peak energy, we 
have inferred that the symmetry energy effect on the isoscalar monopole is mainly due to 
variations of the "Slope"-L around saturation. This is an important result, since it could 
open the possibility of a direct access to this poorly known parameter from monopole data 
changing the neutron number in a fixed isotope. Moreover we will be able to trace back the 
covariant structure of the effective interaction in the isovector channel, of relevance also for 
relativistic heavy ion collisions 22]. 



The isovector dipole response is directly linked to the isospin dependent part of the nu- 
clear EoS. The new result shown here is that the IVGDR energy is decreasing when the 
5 coupling is included, keeping fixed the symmetry energy at saturation. The effect can be 
easily understood from the form of isovector dispersion relations in nuclear matter in a rela- 
tivistic linear response approach. We can then have different IVGDR energies for effective 
interactions that show the same 04 value, but with a different covariant field structure in 
the isovector channel. We note the similarity to the isoscalar case, just exchanging the roles 
of compressibility and isoscalar couplings with symmetry energy and isovector couplings. 

Finally we do not see large effects from the Density Dependence of the effective meson 
couplings. This is mainly due to the smooth behavior around po- Moreover we have shown 
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that some compensation is also coming from dynamical contributions of the rearrangement 
terms. 
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